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ABSTRACT 

Context. Numerical models of the wind-blown bubble of massive stars usually only account for the wind of a single 
star. However, since massive stars are usually formed in clusters, it would be more realistic to follow the evolution of a 
bubble created by several stars. 

Aims. We develope a two-dimensional (2D) model of the circumstellar bubble created by two massive stars, a 40 Mq 
star and a 25 M© star, and follow its evolution. The stars are separated by approximately 16 pc and surrounded by a 
cold medium with a density of 20 particles per cm"^. 

Methods. We use the MPI-AMRVAC hydrodynamics code to solve the conservation equations of hydrodynamics on a 
2D cylindrical grid using time-dependent models for the wind parameters of the two stars. At the end of the stellar 
evolution (4.5 and 7.0 million years for the 40 and 25 M0 stars, respectively), we simulate the supernova explosion of 
each star. 

Results. Each star initially creates its own bubble. However, as the bubbles expand they merge, creating a combined, 
aspherical bubble. The combined bubble evolves over time, influenced by the stellar winds and supernova explosions. 
Conclusions. The evolution of a wind-blown bubble created by two stars deviates from that of the bubbles around single 
stars. In particular, once one of the stars has exploded, the bubble is too large for the wind of the remaining star to 
maintain and the outer shell starts to disintegrate. The lack of thermal pressure inside the bubble also changes the 
behavior of circumstellar features close to the remaining star. The supernovae are contained inside the bubble, which 
reflects part of the energy back into the circumstellar medium. 

Key words. Hydrodynamics - Stars: circumstellar matter - Stars: winds: outflows - Stars: Supernovae: general - ISM: 
bubbles - ISM: supernova remnants 



1. Introduction 

It is well-known that the winds of massive stars cre- 
ate large-scale (> 10 pc) bubbles around their progeni- 
tor s. An analytical mod el of these bubbles was presented 
by I Weaver et al.l (jlQTTT ). Numerical models, which incor- 
porated the time-dependent nature of stellar winds were 
developed b y, for e x ample, iGarci a-Scgura et al. (1 9961 



found in the Rosetta nebula around clu ster NGC 2244 (e.g. 
iPhelDS fc Ladalll997l: EFfc Smithll200,'5h . 

Because of the numerical problems associated with mod- 
eling the bubbles of multiple stars, attempts to account for 
the influence of multiple stars on the evolution of the bub- 
ble have been limite d to the details of the colliding winds 
of bi n ary stars (e.g. Stevens et al.l 19921: Folini Sz Walded 



iPrever et all (I2003D. Ivan Marie et al.l (|2005D . iDwarkadaa 
(|2007l ) and iToala fc Arthuil (|201li rLike the analytical 
model, these models only considered the evolution of a bub- 
ble blown by a single star. However, since massive stars are 
typically formed in star clusters, it is unlikely that any sin- 
gle star could form a bubble on a scale of tens of parsecs, 
without encountering other massive stars. The expanding 
bubble should instead collide with other, similar bubbles 
and m erge as in the case of WN8-WN9h (" Mauerhan erall 
I2OIOI ). A good example of the resulting bubble can be 
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I2OOOI: iPittard fc ParkinI I2OIOI: Ivan Marie et al.l 12011'). or 
dealt with only part of the stellar evolution (Stcgemanj 
'20031). IVelazauez et al.l (|2003l) simulated the collision be- 
tween a supernova and an expanding wind bubble, but with 
a simplified model that assumed the supernova progenitor 
had not previously created a circumstellar bubble through 
its own wind. 

In this paper, we present the result of a two-dimensional 
(2D) numerical simulation that models the evolution of a 
bubble, blown by the winds of two massive stars and com- 
pare these to the bubbles blown by each star individually. 
We follow the evolution of a 25 M0 star that passes through 
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Table 1. Wind and supernova parameters, with tend the 
time from the start of the simulation to the end of each 
evolutionary phase, m the wind mass loss rate, the wind 
veloicty, and M and E the total mass and energy injected 
into the medium. 
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the main sequence (MS) and the red supergiant (RSG) 
phase, and a 40 star that passes through the MS, RSG, 
and Wolf-Rayet (WR) phases. Since the stars are approx- 
imately 16 pc apart, they initially form individual, spher- 
ical bubbles that expand into the surrounding interstellar 
medium (ISM). As the bubbles expand, they collide and 
merge, forming an aspherical bubble, fed by the winds from 
both stars. Over time, the wind parameters change, alter- 
ing the morphology of the bubble, until the stars reach the 
end of their evolution and explode as supernovae. 

We provide movies of the evolution of the bubbles as 
online material (Appendix A). 

2. Numerical setup 

We use th e MP I-AMRVAC hydrodynamics code 
(jKeppens et al.l I2012I ). which solves the conservation 
equations for mass, momentum, and energy on an adaptive 
mesh grid. The effect of r adiative cooling is included 
(Ivan Marie fc Keopensl l2nilh . using a cooling curve for 
solar-metallicity gas generated with the CLOUDY code 
(Ferland et al. 1998; Wang Ye, private communication). 
We set a minimum temperature of 100 K throughout the 
simulation, which limits the amount of compression due to 
radiative cooling, preventing numerical problems. Here we 
neglect the effect of photoionization, which would increase 
the gas temperature within the Stromgren radius to about 
10,000K. 

We set up our simulation as a 2D cylindrically sym- 
metric grid in the r, z-plane. The grid has a basic reso- 
lution of 120x240 gridpoints over 100x200pc. The adap- 
tive mesh has four layers, allowing for an effective resolu- 
tion of 960x1920 gridpoints. At the start of the simula- 
tion, this grid is filled with gas with a constant density of 
10~^^'^gcm~^ and solar metallicity. The stars are defined 
as small spheres {R = 1.5 pc), placed 16.2 pc apart that 
are filled, at the beginning of each timestep, with material 
according to the wind parameters given in Table [TJ The 
duration of the evolutionary phases and the mass-loss rates 
are based on the parameters f o r a 25 M© s tar and a 40 M(7^ 
star u sed by Ivan Marie et al.l (|2004l) and Ivan Marie et al.l 
(|2005D . respectively. W ind velocities reflect those typically 
found in observations (jLamers fc Cassinellilll999l and ref- 



Fig. 1. Logarithm of the density of the circumstellar 
medium in g/cm^ at the end of the evolution of the 40 M0 
star (left) and the 25 M© star (right), if the two bubbles 
never merge. 



erences therein). For the 40 M© star, the mass-loss rate in 
the final phase (the Wolf-Rayet ph ase) has been reduc ed to 
reflect the values obtained by Vink fc de Koteil (|2005[ ). 

When a star reaches the end of its evolution, we fill a 
sphere {R = 0.75 pc) with a constant density and thermal 
energy according to the supernova values given in Table [TJ 
The masses reflect the fin al masses of the stars as well the as 
values calculated by, e.g.. lEldridge fc To"utl (|2004[ ). For sim- 
plicity, we use a supernova energy o f 10^^ erg, which is typ- 
ical of core-collapse supernovae (e.g. Ivan Marie et al.ll2010L 
and references therein) 



3. Results 

3.1. Single stars 

For comparison, we run the evolution of the circumstellar 
bubble for each star individually. Fig. |T] shows the density 
of the circumstellar medium at the end of the evolution 
of each star. The 40 M© star has created a large bubble 
(R'^ 25 pc) that has s wept up a thin, shell containing linear 
thin-shell instabilities (|Vishniac ,1983'). The star has already 
passed through the RSG phase to become a WR star, and 
the W R wind has swept up the remnant s of the earlie r 
phases (jGarcia-Segura et~a l. 1996a; van Ma rie et al.ll2005D . 
The 25 Mq star ends its evolution as a RSG and forms a 
thin shell at the wind termination shock ( ^ 3 pc) . The 
total bubble is smaller (R'^ 23 pc) than for the 40 Mq star, 
though the difference is small because the longer life of 
the smaller star partially compensates for the lower ram 
pressure of its wind. The swept-up shell at the outer edge 
is much thicker. This is the result of the lower ram pressure 
of the wind, which causes the bubble to expand at a lower 
rate, giving the swept-up shell more time to expand under 
its own internal pressure. Because the shell is thicker, it is 
not subject to thin-shell instabilities and only shows the 
beginning of Rayleigh- Taylor instabilities at its inner edge. 
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Fig. 2. Logarithm of the density of the circunistellar medium in g/cm^ after 1.17 and 2.35 Myr (left panel), after 4.42 
and 4.58 Myr (center panel), and after 4.774 and 4.783 Myr (right panel). On the left, the shells of the two bubbles 
encounter each other (left panel, left side) and are destroyed in the process. The remnants are pushed toward the lightest 
of the stars (left panel, right side). On the right, the 40 Mq star (at Y = —2.5 x 20^^cm) goes through the RSG phase, 
which forms a thin shell at the wind termination shock (center panel, left side) and then the WR phase. The WR wind 
sweeps up the RSG wind and the RSG shell is destroyed in the collision. The debris of the collision moves outward into 
the hot shocked gas (central panel, right side). Finally, the 40 Mq star explodes as a supernova, which expands into the 
circumstellar bubble (right panel, left). Once it has hit the outer shell, the supernova bounces back and starts flowing 
inward (right panel, right). 
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Fig. 3. Logarithm of the density of the circumstellar medium in g/cm'^ after 4.86 and 5.60 Myr (left panel), after 6.58 and 
7.0 Myr (center panel), and after 7.01 and 7.07 Myr (right panel). On the left, the mass, ejected by the first supernova 
has spread throughout the bubble, which, due to the injection of energy, expands rapidly. The density inside the bubble 
decreases and the Ray leigh- Taylor instabilities in the outer shell become more pronounced. In the central panel, the 
25 Mq star has reached the RSG phase and forms a shell at the wind termination shock that moves away from the star 
as the thermal pressure in the bubble decreases over time. Finally, on the right, the 25 Mq star explodes as a supernova, 
which sweeps up the RSG shell and fills the bubble. 



3.2. Two stars combined 

Each star initially creates its own circums tellar bubble. 
These bubbles follow the pattern predicted bv lWeaver et al.l 
(|1977f ): The wind expands freely until it hits the termina- 
tion shock. The kinetic energy of the wind is then turned 
into thermal energy creating a volume filled with very hot 
(T ~ 10^ - 10*K) shocked wind material. The high ther- 
mal pressure of the shocked wind material causes it to ex- 
pand, sweeping up a shell of shocked ISM, which moves 
outward. Eventually, the swept-up shells of the two bub- 
bles encounter each other and both are fractured during 
the collision (left panel of Fig. [2]). The remnants are pushed 
into the bubble with the lowest thermal pressure (resulting 
from the weakest wind, and therefore, typically, the lightest 



star) . Once the shells have been fractured, the two bubbles 
merge, creating a single, aspherical bubble in which the 
thermal pressure is equalized as the shocked gas of the two 
stellar winds mixes. The th in outer shell i s unstable and 
shows both linear thin-shell (jVishniadll983[ ) and Rayleigh- 
Taylor instabilities. 

The 40 M0 star, which has the shortest lifespan, leaves 
the main sequence and goes through the RSG and WR 
phases (central panel of Fig. [2]). The RSG wind forms a 
thin shell at the termination shock, which is destroyed by 
the WR wind. Finally, the star explodes as a supernova. 
At this time, the bubble is still quite aspherical, despite 
the outer shell being driven by the thermal pressure in the 
shocked wind. It has a maximum length (along the axis 
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connecting the stars) of 57 pc. At the position of the 40 Mq 
star, it is ~52 pc wide, whereas at the position of the 25 M© 
star the width is only 39 pc. The shape of the outer shell 
still indicates that it is the merger of two spherical bubbles. 

The supernova quickly expands into the low density 
bubble that has been created by the stellar winds, but stops 
when it reaches the contact discontinuity separating the 
shocked wind from the shocked ISM (right panel of Fig. [5]) . 
The injection of energy from the supernova (approximately 
three times the total energy injected so far by the stel- 
lar winds) causes the bubble to expand more rapidly (left 
panel of Fig. [3]). Because the supernova has injected a lot 
of energy, but (relatively) little mass, the density inside 
the bubble decreases quickly during the expansion and the 
Ray leigh- Taylor instabilities in the outer shell grow quickly. 
The bubble is left with only one source of energy, the wind 
from the 25 M© star, which is insufficient to maintain the 
expansion of the bubble. As a result, the thermal pressure 
inside the wind bubble decreases leading to a loss of com- 
pression at the outer shell, which starts to expand because 
of its own internal pressure (left panel of Fig. During 
this phase, the asymmetrical shape of the bubble starts 
to disappear as the isotropic pressure in the shocked wind 
bubble pushes the shell outward equally in all directions. 

Eventually, the 25 M© star reaches the RSG phase, 
which further reduces the energy being added to the bubble. 
A thin shell forms at the wind termination shock (central 
panel of Fig. [3]). This initially happens close to the star 
because of the low ram pressure of the slow RSG wind. 
However, owing to the relatively low pressure in the bub- 
ble, this shell can move farther away from the star than 
would have been possible if the 25 Mq star had been alone 
(i? ~ 6.5 pc, rather than R 4.8 pc when comparing the 
right side of the center panel of Fig. |3] with Fig [ij . 

Finally, the 25 Mq star reaches the end of its evolution 
and explodes (right panel of Fig. [3]), sweeping up the sur- 
rounding RSG wind as well as the RSG shell. The bubble 
has lost most of its asymmetry; it has a length of '^105pc 
and a maximum width of ~104pc. The outer edge of the 
shell has an almost spherical shape. As in the case of the 
first star, the supernova quickly expands into the low den- 
sity shocked wind bubble until it reaches the outer shell, 
then falls back inward. From this moment on, the bubble 
has no energy source. While energy continues to be lost 
from the bubble by thermal conduction, turbulent mixing, 
and radiative cooling, the outer shell starts to push inward 
owing to is own internal pressure, which is considerable be- 
cause of its relatively high density. Assuming a typical speed 
for gas in the shell of 1...2 kms""'^ (the sound speed for gas 
at a 10^ — 10'^ K), the gas in the outer shell would take 
a minimum of 2.5 x 10^ years to reach the center, though 
in practice it would take longer owing to the continuing 
resistance of the hot bubble. 



4. Discussion and conclusions 

The presence of multiple stars inside the circumstellar bub- 
ble influences the evolution of the bubble in several ways. 

The most obvious is in terms of the outer shape of the 
bubble, which is strongly aspherical during the early phase 
of its evolution. This effect is clearly visible until after the 
supernova explosion of the first star. In addition, remnants 
of the collision between the two outer shells of the individual 



bubbles remain visible within the shocked wind bubble until 
after the supernova explosion of the 40 Mq star. 

The collapse of the outer shell, which is caused by a 
decrease in the internal pressure of the bubble, is unique to 
the multi-star scenario. It can only occur when the bubble 
persists for an extended period of time 10^ — 10^ yr) with 
a significantly diminished energy source (the wind of the 
remaining star). The effective result is that the outer shell 
resembles the single-star shell of the star with a stronger 
wind. As long as the 40 Mq star exists, the shell looks like 
a shell created by a single 40 Mq star. Afterwards, it starts 
to resemble the circumstellar shell of a 25 Mq star as it 
becomes relatively thick. However, it is far more instable as 
a result of the instabilities created in the initial evolution. 
This lack of internal pressure also affects the medium close 
to the 25 Mq star: the RSG shell can move further away 
from its progenitor than would have been the case for a 
single-star bubble. 

The supernovae are contained by the outer shell and do 
not break out of the bubble. This is surprising, considering 
the amount of energy in the supernova. However, the mass 
of the outer shell is very high. A bubble with i? = 30pc 
and an ambient medium density of 10~^^'^gcm~^ has a 
5.3x10'* Mq outer shell that forms a very effective barrier 
against a supernova, which is high in energy, but low in 
mass. 

We neglected the effect of photoionization, which would 
give the gas a minimum temperature of about 10,000 K 
during the main sequence and WR phases. During these 
phases, the shells would become thicker owing to their 
higher thermal pressure. This increase in pressure within 
the outer shell would also hasten its expansion after the 
first supernova. 

We conclude that simulating the bubbles around more 
than one star improves our understanding of the evolution 
of the circumstellar bubbles of massive stars. 
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Appendix A: Animations of the evolution of 
circumstellar bubbles 
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Fig. A.l. The first movie (starl) shows the evolution of the bubble created by the 40 star on its own. It follows 
the cvohition from the main sequence, to the RSG stage, and then onto the end of the WR phase. The interval between 
frames is 39112 years. 

Fig. A. 2. The second movie (star2) follows the evolution of the bubble created by the 25 Mq star on its own. It follows 
the evolution to the end of the RSG phase. The interval between frames is 391 years. 

Fig. A. 3. The third movie (wwsn_prcsnl) shows the evolution of the bubble prior to the first supernova. The 40 star 
evolves through the main sequence, the red supergiant phase, and then the Wolf-Rayet phase. The 25 Mq star remains 
on the main sequence. As the two bubble merge, the shells collide and break up. The remnants are pushed toward the 
star with the lowest mass, because it has the weakest wind. The interval between frames is 39 112 years. 

Fig. A. 4. The fourth movie (wwsn_snl) shows the expansion of the first supernova. The 40 star explodes as a 
supernova, which expands within the wind-blown bubble of the two stars. After colliding with the outer shell, the 
supernova remnant falls back toward the center of the bubble. The interval between frames is 391 years. 

Fig. A. 5. The fifth movie (wwsn.postnl) shows the evolution of the bubble between the two supernovae. The second 
star evolves through the remaining main sequence phase and the red supergiant phase. The bubble, which no longer 
derives energy from the two stars, loses the thermal pressure necessary to sweep up the outer shell, which starts to fall 
apart. The interval between frames is 39 112 years. 

Fig. A. 6. The sixth movie (wwsn_sn2) shows the expansion of the second supernova. The second star explodes as a 
supernova, which expands inside the bubble. The interval between frames is 391 years. 



